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Abstract We present the results of velocity measurements obtained by ultrasonic Doppler velocimetry and local 
potential probes in the flow of GalnSn eutectic melt driven by a two phase inductor in a cylindrical container. This 
type of flow is expected in a recent modification to the floating zone technique for the growth of small-diameter 
single intermetallic compound crystals. We show that the flow structure can be changed from the typical two toroidal 
vortices to a single vortex by increasing the phase shift between the currents in the two coils from to 90 degrees. The 
latter configuration is thought to be favourable for the growth of single crystals. The flow is also computed numerically 
and a reasonable agreement with the experimental results is found. The obtained results may be useful for the design 
of combined two-phase electromagnetic stirrers and induction heaters for metal or semiconductor melts. 

1 Introduction 

The present work is concerned with a physical modelling of the melt flow in a recent modification to the floating zone 

technique for the growth of small-diameter single intermetallic compound crystals [3]. The quality of crystals grown 

by this method depends on the growth conditions, particularly on the shape of the solidification front. The growth 

of single crystals usually requires a convex solidification front pQ. The shape of this front can strongly be affected by 

the convective heat transport in the melt. A standard single-phase induction heater usually drives a radially inward 

jet of hot melt at the middle of the floating zone [2j. At the centre of the floating zone, the radial jet splits into two 
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nearly symmetric axial jets, which further stream towards both solid-liquid interfaces. The hot axial jet impinging at 
the centre of the growth interface renders it concave so promoting a polycrystalline growth. To overcome this adverse 
effect, we recently proposed a two-phase inductor, which is able to reverse the direction of the electromagnetically 
driven melt flow so that the solidification front becomes mostly convex [3]. 

The inductor consists of two coils that are fed by alternating currents of the same amplitude but a 90 degrees 
phase shift between them. In the original design for radio-frequency induction heating, the phase lag is achieved by 
short circuiting the secondary coil through a variable capacitor that is tuned so that the resonance frequency of the 
secondary circuit coincides with the frequency of the power supply [5]. In this way a component of the magnetic field 
is created that travels towards the coil with the phase lag so dragging the liquid metal along it. As a result, the outer 
layers of the hot melt are driven towards the periphery of the solidification front while the cooled down melt returns 
in the centre [6]. A similar multiphase stirrer for low conductivity melts has been developed by Ernst et al. [7], who 
also used it for refining the grain size in gold alloys |5j. 

Until now there have been no direct experimental observations of the effect of the phase shift in a two-coil inductor 
on the structure of ensuing fluid flow. In this paper, we present comprehensive results of velocity measurements in 
a model experiment using GalnSn eutectic alloy driven by a simple two-phase stirrer in a cylindrical container. The 
evolution of the flow pattern has been determined depending on the phase shift between the two coils which were 
connected to separate power supplies with externally controllable phase shifts. The experiments were also modelled 
numerically and numerical results compared with the measurements. 

The paper is organized as follows. In the following two sections, experimental setup and measurement techniques 
are described. Experimental results are presented in Sec. [3] and compared with numerical results in Sec. [5] The paper 
is concluded by a summary in Sec. [6] 

2 Experimental setup 

The main element of the experimental setup shown in Fig. Q]is an open top cylindrical container of 50 mm in both the 
diameter and heigh (unity aspect ratio) made of a 2 mm thick Polycarbonate. The container was filled with GalnSn 
alloy at eutectic composition with 67% Ga, 20% In and 12.5% Sn, and kept in a water jacket at 15 °C temperature. 
The relevant physical properties of GalnSn are listed in Table [T] The top surface of the liquid metal was kept open 
for measurements. The Joule heat generated in the liquid metal was efficiently removed by passing the cooling water 
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Figure 1 Experimental setup in cross-section (a) and 3D (b) views showing the inner container of diameter D = 50 mm with 
liquid metal (1), the outer container with cooling water (2), two single-winding coils of diameter D iud = 88mm and vertical 
separation d = 44 mm (3), one of the six cooling water intakes and outlets distributed uniformly along the circumference (4), 
water jacket (5), cooling water partitioning walls (6), water intakes and outlets for coil cooling (7), GalnSn alloy (8), coil 
terminals (9). 

through a labyrinth of thin partitioning walls, which were uniformly distributed along the side wall of the container. 
The container was placed inside an inductor consisting of two circular copper coils of diameter Di n d = 88 mm, axial 
separation d = ^D in ^ and a rectangular 8 x 9 mm 2 cross-section with an inner cooling channel. To minimize the 
perturbation of the magnetic field by the coil terminals, the latter were made of a 2 mm thick copper strip welded into 
each copper winding with a 1 mm gap, which was filled with Teflon. Thus the current in each coil was interrupted only 
over this small gap which constituted less than 0.4% of the total coil circumference. To minimize the superposition 
of perturbations, the terminals of each coil were placed across the diameter of the container as shown in Fig. Hp)). 

The coils were powered as follows. First, two source signals with preset amplitudes, frequencies and relative phases 
were generated digitally by a computer program and then converted into analogue form by a D/A converter. Second, 
the analogue signals were amplified by a two-channel power amplifier. Two toroidal transformers were used to match 
the output impedance of the amplifier with that of the inductor. This allowed us to reach the RMS current up to 
600 A in each coil corresponding to a maximum magnetic field of 6.3 mT at the centre of the container. The AC 
frequency was set to / = 400 Hz which is close to the optimal one that produces maximal flow velocity at a fixed 
current amplitude in this setup [T4] . 
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density (kg/m 3 ) 


6360 


electrical conductivity (S/m) 


3.2 x 10 6 


kinematic viscosity (m 2 /s) 


3.4 x icr 7 


container diameter X height (m) 


0.05 x 0.05 


diameter of coils (m) 


0.088 


distance between coils (m) 


0.044 


current amplitude (A) 


100 .. . 600 


AC frequency (Hz) 


400 



Table 1 Physical properties of the GalnSn alloy and experimental parameters. 

The velocity distribution in the liquid metal was measured by two methods: pulsed ultrasound Doppler velocimetry 
(UDV) and potential probe with local magnetic field, which both are briefly described below. 



3 Measurement techniques 

3.1 Ultrasound Doppler velocimetry 

Ultrasound Doppler pulsed velocimetry (UDV) is a widely employed technique for measuring velocities in liquid 
metal flows [9 ll01lllj . The method uses the phase shift of the ultrasonic wave scattered by microscopic particles in 
the liquid to determine the distance to the particles. By analysing the time correlation of the signals acquired at 
different instants of time the device reconstructs the velocity of particles along the ultrasonic beam. In this study 
we employed DOP2000 Velocimeter model 2125 ('Signal Processing SA', Switzerland) to measure the axial velocity 
profiles along the height of the container in the flow of GalnSn driven by a two-phase AC magnetic field. Measuring 
velocities of a few centimetres per second over a relatively small height of the container required the highest available 
ultrasound frequency of 10 MHz. 

The ultrasound wave was introduced into the liquid metal directly through its free surface. A special care was taken 
to ensure a good wetting of the UDV transducer by the liquid metal which was essential for the acoustic coupling. 
The vertical velocity profiles were measured at the following 9 radial positions: r = 0, ±5, ±10, ±15, ±21 mm. The 
radial size of the UDV probe precluded measurements at larger r, i.e. closer to the side wall. At each radial position 
300 samples of axial velocity profiles were recorded during a 35 s interval with ~ 8.5 Hz sampling rate and a 0.3 mm 
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Figure 2 Potential probe consisting of a small permanent magnet with two pairs of electrodes. The horizontal and vertical 
pairs of electrodes measure the opposite, that is the vertical and horizontal, respectively, velocity components. 

spatial resolution in the axial direction. These data were used to find both the time average of the velocity and its 
standard deviation, i.e. the intensity of velocity pulsations, at each of 9 x 170 spatial measurement points over the 
vertical cross-section of the container. Since the measurements at different r were taken at different times, it was not 
possible to follow the temporal evolution of the flow over the whole cross-section. 

3.2 Potential probe 

Potential probe, also known as the Vives probe [12, 13 , is another commonly used technique for local velocity meas- 
urements in low-melting-point liquid metals. The probe shown in Fig. [2] consists of a small permanent magnet with 
two pairs of small electrodes. The liquid metal flow passing between the electrode pairs in the magnetic field of 
the probe induces a difference of the electric potential between the electrodes proportional to the velocity in the 
respective direction. Besides the low cost, the main advantage of the potential probe over a single-channel UDV is 
the possibility to measure simultaneously two components of local velocity with a high temporal sampling rate. How- 
ever, there are also several disadvantages. First, very low voltages, typically about 0.1 . . . 0.2 /iV/(mm/s), need to 
be measured. Second, quite long measurement times are required even at moderate spatial resolutions to acquire the 
velocity distribution over the whole cross-section of the container. Third, the temperature of the liquid metal has to 
be kept sufficiently uniform to avoid thermo-electrically induced potential perturbations. In the present experiment, 
the potential probe was used only for a few measurements to check the UDV results. 
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Figure 3 Calibration curve for the potential probe with horizontal bars showing the RMS of fluctuations. 



The probe used in this experiment to measure the axial and radial velocity components was made of a cylindrical 
NdFeB permanent magnet of 2 mm in both diameter and height, which was hold by a stainless steel tube of 1 mm in 
diameter. The probe was calibrated using an annular channel filled with GalnSn which was rotated about its centre 
with a well defined angular velocity. Since the sensitivity of the probe is non- linear at low velocities (see Fig. [3}, 
calibration was carried out for the whole range of the expected velocities. The calibration data were fitted with two 
mathematical expressions: a 'sigmoid law' for low velocities and a linear one for higher velocities. These expressions 
were further used in the acquisition software to convert the voltage measurements into velocities. 

Measurements were taken at 121 locations in the vertical cross-section of the container. At each location the 
axial and radial velocity components were recorded during a period of 100 seconds with a 50 Hz sampling rate using 
high-resolution 24-bit 'Data Translation' DT9821 acquisition board. The internal pre-amplifier of the board with the 
gain of 64 ensured a 20 nV 'noise-free' resolution in voltage at about the same RMS value of the noise. In order to 
prevent a thermally induced zero drift, the board was kept in a thermo-box. 



4 Experimental results 



4.1 UDV measurements 



The basic effect of the phase shift between the currents in the top and bottom coils on the flow structure is illustrated 
by Fig. U which shows the isolines of the vertical velocity component for the magnetic field travelling upwards (a) 
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Figure 4 Isolines of the axial velocity component for the magnetic field travelling downwards (a) and upwards (b) which 
correspond to the current in the upper coil lagging and leading, respectively, in phase by 90° relative to the current in the 
bottom coil. In both coils the effective current is the same and equal to / = 350 A. 

no-i 

■ Experimental data (Kith StDev vol.) 
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Figure 5 Measured by the UDV and computed maximal absolute value of the axial velocity component versus the effective 
current / at the phase shift fixed to Aip = 90° . 

and downwards (b). Note that the velocities are recorded as positive when the flow is directed away that is downward 
from the UDV probe which is placed at the top of the liquid metal. Also note that the velocity could be measured 
by the UDV only up to |r| < 21 mm and is then extrapolated to the region near the side wall where it is supposed to 
obey the no-slip condition. In both cases the currents in the top and bottom coils alternate with the same frequency 
of / = 400 Hz and have the same effective values equal to I = 350 A. The only difference is that in the first case 
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the current in the upper coil lags in phase by 90° while in the second case it leads by the same amount relative to 
the current in the bottom coil. As the lagging current follows the leading one, the magnetic field effectively travels 
from the former to the latter. And as it does so, it drags the liquid metal near the side wall with it. Respectively, 
the magnetic field travelling upward gives rise to regions of negative axial velocity at the side walls, while there 
is a larger region of positive velocity in the central part of the container where a downward axial jet forms. The 
highest axial velocity is attained in the jet at the symmetry axis (r = 0) somewhat downstream from the mid-height 
of the container. Switching the sign of the phase shift and thus the vertical direction in which the magnetic field 
travels swaps both the direction of the flow and its pattern about the mid-height of the container. The slight axial 
asymmetry of the flow pattern observable in both cases may be due to the magnetic field perturbations by the current 
connections. 

Two series of measurements were performed using the UDV. In the first series, the effective current I in both coils 
was increased simultaneously from 100 to 600 A while the phase shift between them was kept fixed at Aip = 90°. 
The maximal axial velocity, which for this case is plotted in Fig. [5] is seen to increase with the current almost linearly 
from 15mm/s to 92 mm/s. Along with the velocity grow also its RMS fluctuations which characterize the level of 
turbulence in the flow. The flown pattern remains basically similar to that shown in Fig. [4] (left panel) with a slight 
axial asymmetry. 

In the second series, the effective currents in both coils were kept fixed at I = 350 A while the phase shift Atp 
between them was increased from 0° to 180°. Figure [S] shows the isolines of the axial velocity component measured 
at four different phase shifts Aip in addition to that of Ap = 90°, which is plotted in Fig. HJa). The first case with 
Aip = 0° shown in Fig.|BJa) corresponds to a standard single-phase AC magnetic field. This field gives rise to a nearly 
compressing electromagnetic force with a maximum at the mid-height between the coils, where it drives the liquid 
metal radially inwards. As a result the typical two- vortex structure, which can be inferred from Fig. [HI a), arises. As 
seen in Fig. [6jb) , a phase shift of Aip = 30° is sufficient to break this two- vortex balance in favour of the bottom one 
which is enhanced by the arising upward axial force while the top vortex is suppressed by it. The axial flow keeps 
increasing with the phase shift up to Aip = 90° where according to Fig. [7] it attains a maximum and then starts to 
decrease when the phase shift is increased further. The lowest axial velocity is attained at Aip = 180° which again 
corresponds to a single-phase AC magnetic field. In contrast to the previous case, the currents in both coils now 
flowing in the opposite directions tend to cancel out each other's magnetic field. Thus, the electromagnetic pinching 
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Figure 6 Isolines of the axial velocity component for various phase shifts Atp between the currents in top and bottom coils 
(a) 0°, (b) 30°, (c) 120°, (d) 180°. The case of Atp = 90° can be seen in Fig. Ha) 

force cancels out in the mid-plane between the coils where the liquid is now pushed by pressure gradient radially 
outwards. This gives rise to a two- vortex structure shown in Fig. H^d), which is similar to that for Atp = 0° seen in 
Fig. (6{ a) except for the opposite flow direction and a twice lower velocity. 

4.2 Potential probe measurements 

Potential probe was used only in one case for / = 350 A and Atp = 90°, which was then compared with the UDV 
measurements and numerical results. The measured velocity field and the isolines of axial velocity component are 
shown in Fig. [8] The latter is seen to agree well with the corresponding UDV measurements shown in Fig. |4ja). 
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Figure 7 Axial velocity profiles along the symmetry axis (r = 0) at different phase shifts between the currents in both coils 
(a). The highest axial velocity along the symmetry axis and the associated RMS velocity fluctuations versus the phase shift 
between the currents (b). 
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Figure 8 The velocity field in the vertical cross-section (a) and the isolines of the axial velocity component (b) measured by 
the potential probe at / = 350 A and A(p = 90° . 



Even the slight axial asymmetry is reproduced. This probe could reach the radial positions up to r = ±23 mm so 
approaching the side wall closer than the UDV well into the region of ascending flow. The full flow pattern shown 
in Fig. [5] includes an extrapolation of the velocity distribution up to the side walls with no-slip boundary condition 
applied. As seen in Fig. [9] there is also a good agreement between the maximal values of axial velocities measured 
by the potential and UDV probes. 
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5 Comparison with numerical results 



The problem was also modelled numerically using the finite element code ANSYS in axially symmetric approximation. 
First, the ANSYS Emag solver was used to solve the electromagnetic part of the problem. At this stage, the effect of 
liquid metal flow on the AC magnetic field was neglected following the inductionless approximation which is commonly 
used in liquid metal magnetohydrodynamics. The cross-section of metal was discretizied with 5000 mapped elements. 
The distribution of the magnetic field was computed in both the liquid metal and the surrounding areas for various 
coil currents and their phase shifts. From those data the distribution of the electromagnetic force in the liquid metal 
was computed and then passed on as a body force to the FLOTRAN hydrodynamic solver of ANSYS. Using this 
solver with a standard k-e model of turbulence and no-slip boundary conditions at the container walls the flow of 
liquid metal driven by AC magnetic field was computed. 

First, as seen in Fig. [5] the numerical results demonstrate a good agreement for the maximal values of the axial 
velocity measured by the UDV technique at various currents and a phase shift fixed to Aip = 90°. There is also 
a reasonable agreement between the measured and computed profiles of the axial velocity along the height of the 
container at various radii, which are plotted for I = 350 A and A<p = 90° in Fig. 1101 The slight axial asymmetry in 
the measured velocity distribution, which may be due to the current connectors, was eliminated in the comparison 
by taking the mean value between positive and negative radii. 
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Espenmenlal data: 
n r=0 
Q 1 = 5 rrrn 
a 1 = 10 rrm 
v i = 15 rrm 
<> i = 21 rrm 

Numerical data: 

r = Drrm 

i = 10 mn 

r = 15 run 
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Figure 10 Measured and computed profiles of the axial velocity along the height of the container at various radii for I = 350 A 
and Aip = 90° . 

6 Conclusion 

The velocity distribution in the flow of GalnSn eutectic alloy driven by a two phase inductor in a cylindrical container 
was measured by ultrasonic Doppler velocimetry and local potential probes. The flow was also computed numerically 
and a reasonable agreement with the experimental results was found. It was demonstrated that the flow structure 
can be modified from the typical two toroidal vortices into a single vortex when the currents in the two coils are 
made to alternate with a certain phase shift. The vortex is strongest for a given current amplitude when the phase 
shift between the currents in both coils is about 90°. The obtained results may be useful for the design of combined 
two-phase electromagnetic stirrers and induction heaters for metal or semiconductor melts in applications where a 
well-defined stirring is required. 

References 

1. Rudolph P (2010) Defect formation during crystal growth from the melt. In: Dhanaraj G, Byrappa K, Prasad V, Dudley 
M. (eds) Springer Handbook of Crystal Growth, 1st edn, Springer, Berlin, pp 159—201. 

2. Hermann R, Priede J, Behr G, Gerbeth G, Schultz L (2001) Influence of growth parameters and melt convection on the 
solid-liquid interface during RF-floating zone crystal growth of intermetallic compounds. J Cryst Growth 223:577—587 

3. Priede J, Gerbeth G, Hermann R, Behr G, Schultz L and Uhlemann H-J (2007) German Patent DE10328859B4 

4. Hermann R, Behr G, Gerbeth G, Priede J, Uhlemann H-J, Fischer F, Schultz L (2005) Magnetic field controlled FZ single 
crystal growth of intermetallic compounds. J Cryst Growth 275:1533-1538 




Velocity measurements in the liquid metal flow driven by a two-phase inductor 



13 



5. Priede J, Gerbeth G (2012) Analysis of two-phase RF inductor for the floating-zone crystal growth, Electr Eng (submitted) 

6. Rudolph P (2008) Travelling magnetic fields applied to bulk crystal growth from the melt: The step from basic research to 
industrial scale. J Cryst Growth 310: 1298-1306 

7. Ernst R, Perrier D, Brun P, Lacombe J (2005) Multiphase electromagnetic stirring of low conducting liquids. COMPEL 
24:334-343 

8. Ernst R, Hamburger J, Gamier C, Mangelinck-Noel N, Ramoni P (2005) Grain size reduction by electromagnetic stirring 
inside gold alloys. Eur Phys J Appl Phys 30:215-222 

9. Takeda Y (1987) Measurement of velocity profile of mercury flow by ultrasound Doppler shift method. Nucl Techn 79: 
120-124 

10. Takeda Y, Kikura H (2002) Flow mapping of the mercury flow, Exp. Fluids 32: 161-169. 

11. Cramer A, Zhang C, Eckert S (2005) Local flow structures in liquid metals measured by ultrasonic Doppler velocimetry. 
Flow Meas Instrum 15: 145-153 

12. Ricou R, Vives C (1982) Local velocity and mass transfer measurements in molten metals using an incorporated probe. Int. 
J. Heat Mass Transfer 25: 1579-1588 

13. Cramer A, Varshney K, Gundrum Th, Gerbeth G (2006) Experimental study on the sensitivity and accuracy of electric 
potential local flow measurements. Flow Meas Instrum 17: 1-11 

14. Fautrelle Y R (1982) Analytical and numerical aspects of the electromagnetic stirring induced by alternating magnetic 
fields. J Fluid Mech 102: 405-430 



